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Abstract 

The search for topologically non-trivial states of matter has become 
an important goal for condensed matter physics. Recently, a new class 
of topological insulators has been proposed. These topological insula- 
tors have an insulating gap in the bulk, but have topologically protected 
edge states due to the time reversal symmetry. In two dimensions the 
helical edge states give rise to the quantum spin Hall (QSH) effect, in 
the absence of any external magnetic field. Here we review a recent the- 
ory which predicts that the QSH state can be realized in HgTe/CdTe 
semiconductor quantum wells. By varying the thickness of the quan- 
tum well, the band structure changes from a normal to an "inverted" 
type at a critical thickness dc. We present an analytical solution of the 
helical edge states and explicitly demonstrate their topological stability. 
We also review the recent experimental observation of the QSH state in 
HgTe/(Hg,Cd)Te quantum wells. We review both the fabrication of the 
sample and the experimental setup. For thin quantum wells with well 
width dqw < 6.3 nm, the insulating regime shows the conventional be- 
havior of vanishingly small conductance at low temperature. However, for 
thicker quantum wells (dgw > 6.3 nm), the nominally insulating regime 
shows a plateau of residual conductance close to 2e^ /h. The residual con- 
ductance is independent of the sample width, indicating that it is caused 
by edge states. Furthermore, the residual conductance is destroyed by a 
small external magnetic field. The quantum phase transition at the crit- 
ical thickness, dc = 6.3 nm, is also independently determined from the 
occurrence of a magnetic field induced insulator to metal transition. 
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1 Introduction 



Our work on the quantum spin Hall (QSH) effect is motivated both by the 
quest for spin based electronic devices and by the search for topologically non- 
trivial states of matter. Unlike logic devices based purely on charge current, 
quantum spintronic devices integrate information processing and storage units, 
could operate with low power consumption and perform reversible quantum 
computations [1,2]. Recently, the theoretical prediction of the intrinsic spin Hall 
effect [3,4] has generated great interest in the field of spintronics, since this effect 
allows for direct electric manipulation of the spin degrees of freedom without 
a magnetic field, and the resulting spin current can flow without dissipation. 
These properties could lead to promising spintronic devices with low power 
dissipation. The spin Hall effect has been observed recently in both n [5] and p [6] 
doped semiconductors, however, it is still unclear if the underlying mechanism 
is intrinsic [3,4] or extrinsic [7,8]. 

Beyond the potential technological applications, the intrinsic spin Hall effect 
has guided us in the search for new and topologically non-trivial states of matter. 
The quantum Hall (QH) state gives the first, and so far the only example of 
a topologically non-trivial state of matter, where the quantization of the Hall 
conductance is protected by a topological invariant [9]. Given the fundamental 
importance of topological quantization in physics, it is highly desirable to search 
for quantum states of matter characterized by non-trivial topological properties 
similar to, but distinct from the QH state. Soon after the theoretical prediction 
of the intrinsic spin Hall effect in doped semiconductors, a rather dramatic 
prediction was made that the intrinsic spin Hall effect could also be realized 
in insulators [10]. Subsequently, the QSH state was independently proposed in 
graphene [11] and in strained semiconductors [12]. The QSH insulator state is 
invariant under time reversal, has a charge excitation gap in the 2D bulk, but 
has topologically protected gapless edge states that lie inside the bulk insulating 
gap. The edge states have a distinct helical property: two states with opposite 
spin-polarization counter-propagate at a given edge [11, 13, 14]; for this reason, 
they are also called helical edge states. The edge states come in Kramers' 
doublets, and time reversal symmetry ensures the crossing of their energy levels 
at special points in the Brillouin zone. Because of this energy level crossing, 
the spectrum of a QSH insulator cannot be adiabatically deformed into that 
of a topologically trivial insulator without helical edge states; therefore, in this 
precise sense, the QSH insulators represent a topologically distinct new state 
of matter. The topological properties of the QSH state are mathematically 
characterized by a Z2 topological invariant [15]. States with an even number 
of Kramers' pairs of edge states at a given edge are topologically trivial, while 
those with an odd number are topologically non-trivial. 

While the initial proposal of the QSH state in graphene [11] provided an 
interesting theoretical toy model, it was soon shown to be unrealistic since the 
spin-orbit gap in this system is extremely small [16, 17]. Bernevig, Hughes and 
Zhang [18] initiated the search for the QSH state in semiconductors with an "in- 
verted" electronic gap, and predicted a quantum phase transition in HgTe/CdTe 
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quantum wells as a function of the thickness dgw of the quantum well. The 
quantum well system is predicted to be a conventional insulator for dqw < dc 
and a QSH insulator for dgw > dc, with a single pair of helical edge state. In this 
paper, we review the basic theory of the QSH state in the HgTe/CdTe system, 
and provide explicit and pedagogical discussion of the helical edge states. We 
also review recent experimental observations of the QSH effect in this system, 
through detailed discussions of the sample preparation and the experimental 
setup. 

2 General Properties of (Hg,Cd)Te Quantum Wells 

In this section we will review the bandstructure of bulk HgTe and CdTe and 
present a simple model which is useful in describing the physics of the relevant 
subbands of a HgTe/CdTe quantum well. HgTe and CdTe bulk materials have 
the zinc-blende lattice structure. This lattice has the same geometry as the 
diamond lattice, i.e., two interpenetrating FCC lattices shifted along the body 
diagonal, but with a different atom on each sublattice. The presence of two 
different atoms removes the inversion symmetry of the crystal and reduces the 
point group symmetry from Oh (cubic) to Td (tetrahedral) . However, even 
though the inversion symmetry is explicitly broken, this only has a small effect 
on the physics. Except for the later discussions on the magneto-resistance in 
the QSH regime, we will ignore the bulk- inversion asymmetry (BIA). 

For both HgTe and CdTe, the important bands near the typical Fermi level 
are close to the F-point in the Brillouin zone, and they are the s-type band (Fg) 
and the p-type band which is split to a J = 3/2-band (Fg) and a, J = 1/2-band 
(Fy) by spin-orbit coupHng. CdTe, as shown in Fig. [1] (a), has a band order 
similar to GaAs with an s-type (Fg) conduction band and p-type valence bands 
(F8,F7) which are separated from the conduction band by a large energy gap 
of ~ 1.6 eV. HgTe as a bulk material can be regarded as a symmetry-induced 
semi-metal. Its negative energy gap of —300 meV indicates that the Fs band, 
which usually forms the valence band, is above the Fg band. The light-hole 
bulk subband of the Fg band becomes the conduction band, the heavy-hole bulk 
subband becomes the first valence band, and the s-type band (Fg) is pushed 
below the Fermi level to lie between the heavy-hole subband and the spin-orbit 
split off band (Fy). Based on this unusual sequence of the states, such a band 
structure is called inverted. Due to the degeneracy between the heavy-hole and 
light-hole bands at the F point, HgTe is a zero-gap semiconductor. Note that we 
have introduced the cubic group representation labels Fg,F7,F8 for bands near 
the F-point because various orderings of these bands in the two materials are 
cumbersome to keep track of and these symbols explicitly convey the symmetry 
of the bands. We will use these symbols for the remainder of the paper instead 
of referring to s-like, and p-like. 

When HgTe-based quantum well structures are grown, the peculiar proper- 
ties of the well material can be utilized to tune the band structure in a unique 
way. For wide QW layers, the confinement is low and the band structure remains 
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Figure 1: (Upper) Bulk bandstructure for HgTe and CdTe (Lower) Schematic 
picture of quantum-well geometry and lowest subbands for two different thick- 
nesses. 

"inverted" . However, the confinement energy increases, when the well width is 
reduced. Thus, the energy states will be shifted and, eventually, the energy 
bands will be aligned in a "normal" way, if the QW width falls below a critical 
size dc- We can understand this heuristically as follows: for thin quantum- wells 
the heterostructure should behave similarly to CdTe and have a normal band 
ordering, i.e., the bands with primarily Fg symmetry are the conduction sub- 
bands and the Fg bands contribute to the valence subbands. On the other hand, 
as the quantum-well thickness is increased we would expect the material to look 
more and more like HgTe which has its bands inverted. So as the thickness 
increases we expect to reach a critical thickness where the Fg and Fg subbands 
cross and become inverted with the Fg bands becoming conduction subbands 
and Fg becoming valence subbands. This is illustrated in Fig. [1] (b). The shift 
of the energy levels with QW width can be seen in Fig. [2l For the band struc- 
ture, self-consistent Hartree calculations have been performed using a 8 x 8 k-p 
model [19]. The notation of the subbands as heavy-hole (_ff)-like and electron 
(£^)-like is according to the properties of the respective wave functions [20]. The 
light-hole- like subbands are energetically remote {E < —100 meV) so that they 
are not depicted in Fig. [2] The transition from a normal band alignment to an 
inverted one can clearly be seen in this figure. For a thin QW layer the quan- 
tum confinement gives rise to a normal sequence of the subbands, i.e., hole- like 
bands form the valence band and electron states are in the conduction band. 
In contrast, a more complex sequence of the energy states is obtained when the 
QW width dqw exceeds a critical value dc ~ 6.3 nm. In this inverted regime, 
the Hi band lies above the Ei subband of the valence band and, consequently, 
is now the lowest conduction subband. We will now develop a simple model and 
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Figure 2: The energy of the states in the quantum well are shown as a function 
of the width of the HgTe QW layer. 

discuss why, immediately after this crossing, we expect to find a time-reversal 
invariant bulk insulating state with topologically protected edge states, i.e., the 
QSH state. 



3 Effective Model of the QSHE in HgTe/CdTe 
Quantum Wells 



The effective model of semiconductors whose outer bands are s-like and p-\ike 
near the F-point is the 8-band Kane model [21]. This is a model derived from 
k-p perturbation theory for bulk materials. To calculate the subband spec- 
trum we must write down a Hamiltonian which is inhomogeneous along the 
quantum well growth direction. In Ref. [18] we used this method to derive an 
effective 4-band model for the subbands that participate in the inversion cross- 
ing. Rather than repeat this calculation we will give some simple symmetry 
arguments about the form of the model. Under the cubic symmetry group the 
relevant Fg-band transforms antisymmetrically under "cubic" -parity because it 
is an anti-bonding band (in CdTe and HgTe) , even though it has an s-wave like 
character. The Fg band transforms symmetrically since it is a bonding band. 
The "conduction" subbands Ei,E2 . . . and heavy-hole Hi , H2 . . . subbands all 
alternate their parity under z — > — z inversion. For our purposes we neglect the 
bonding combination of the Fg and light-hole bands whose subbands are denoted 
Li, L2 . . . because these bands are well-separated in energy from the bands in 
which we are interested. We will be focussed on only the lowest subbands and as 
mentioned before, the Ei subband actually contains components of the Fg band 
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and the Tg light-hole band i.e. \Ei,±l/2) = a{z)\-f^, ±1 /2) + l3{z)\j^, ±1/2). 
In order to preserve the transformation of the Ei subbands under inversion of 
the growth (z)-direction the envelope function a{z) must be even under z —z 
while the envelope function f]{z) piece must be odd. 

The relevant subbands, Ei and Hi, under our assumption of inversion sym- 
metry, must be doubly degenerate since time-reversal symmetry is present. We 
will label the states with the basis \Ei+),\Ei~),\Hi+),\Hi—), where \Ei±) 
and \Hi±) are the two sets of Kramers' partners. The states \Ei±) and \Hi±) 
have opposite transformations under parity, so a Hamiltonian matrix element 
that connects them must be odd under parity since we assumed our Hamilto- 
nian preserves inversion symmetry. Thus, to lowest order in k, {\Ei+), 
and {\Ei—),\Hi—)) will each be coupled generically via a term linear in k. 
The \Hi+)) heavy-hole state is formed from the spin-orbit coupled p-orbitals 
\Px + ipy,^)), while the \Hi—)) heavy-hole state is formed from the spin-orbit 
coupled p-orbitals | — {px — ipy), i))- Therefore, to preserve rotational symme- 
try around the growth-axis (z-axis) the matrix elements must by proportional 
to k± = kx ± iky. The only terms allowed in the diagonal elements are terms 
that have even powers of k including fc-independent terms. The subbands must 
come in fully degenerate pairs so there can be no matrix elements between the 
-|- state and the — state of the same band. Finally, we note that if non-zero 
matrix elements coupled |£^i-|-), ) or \Ei—),\Hi+) then this would induce 
a higher-order process that would couple the ± states of the same band, thus 
splitting the degeneracy. So these matrix elements are forbidden as well. These 
simple arguments lead us to the following model 



where we have used the basis order {\Ei+) ,\Hi+) , \Ei-), A,B,C,D,M 
are material parameters that depend on the quantum well geometry, and we have 
chosen the zero of energy to be the valence band edge of HgTe at k = as in 
Fig-Ill Ti. is equivalent to two copies of the massive Dirac Hamiltonian but with 
a fc-dependent mass M{k). 

For the purposes of studying the topological properties of this system, as 
well as the edge states, we will work with a lattice regularization of this model 
which will give us the dispersions over the entire Brillouin zone torus, i.e., a 
simplified tight-binding representation. Since all of the interesting physics at 
low energy occurs near the F-point, the behavior of the dispersion at high-energy 
is not important. Thus, we can choose a regularization to make our calculations 



n 




(1) 



h{k) 
e{k) 
da{k) 
M{k) 



e(fc)l2x2+da(fc)fT" 

C-D{kl + kl) 
{Ak,,-Aky,M{k)) 
M-B{kl + kl). 



(2) 
(3) 
(4) 
(5) 
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simple. This simplified lattice model is 



h{k) 
e{k) 
da{k) 
M{k) 



n = 



( Hk) \ 

I h*{-k) ) 

e(fc)l2x2 +rfa(fcK 



(6) 



C — 2Da ^(2 — cos(fca;a) — cos(fcj,a)) 
(^Aa^^ sin{kxa), —Aa~^ sin(fcj^a), M(fc)) 
AI — 2Ba^^ (2 — cos(fca;a) — cos(A:j^a)) . 



It is clear that near the F-point the lattice Hamiltonian reduces to the k-p 
Hamiltonian in Eq. [1] For simplicity, below we take the lattice constant a — I, 
which corresponds to a redefinition of the units. 

The bulk band structure of this model has the energy spectrum 



This spectrum gives a bulk insulator as long as \/dada ^ (and no indirect gap 
is formed) and there is a gapless spectrum at a phase transition critical point 
when this quantity vanishes. The phase transition and topological properties do 
not depend on e(fc) so we will ignore this term for now. The lattice Hamiltonian 
has several critical points, i.e., several points in (k^, ky, M)-space where da = 0, 
but there is only one that occurs in the low-energy spectrum near the F-point. 
Since the model is only valid for small k we will confine ourselves to this single 
critical point. We can see from the form of da that this critical point will occur 
at kx = ky = and M = (assuming B is finite). For M close to 0, the gap 
is minimal around F point, and the system is described by two fiavors of Dirac 
fermions with mass AI and —AI. The Hall conductance of the massive Dirac 
model jumps by a quantum (e^/h) when the mass changes sign [22]. That is, 
aniAI -> 0+) - aniAI 0") = e^/h. For each 2x2 block of the lattice 
model (Eq. (5]) , the Hall conductance can also be calculated explicitly through 
the Kubo formula, which leads to an — cie^ /h, with ci G Z the following 
Pontryagin winding number [23,24]: 



Direct numerical calculations show that ci = for AI/2B < and ci = ±1 
for < AI/2B < 2, where +{~) corresponds to the h{k) {h*{-k)) block. Con- 
sequently, the system with M/2B > is equivalent to two independent quan- 
tum Hall systems with opposite Hall conductance, which guarantees a pair of 
counter propagating edge states on each boundary of the system. (The condi- 
tion AI < AB is always satisfied for the physical system since the gap 2M is 
always small.) Once such a pair of edge states is established on the boundary 
they are robust under time- reversal invariant perturbations, and are protected 



E± 



e{k) ± y/ dad, 



(7) 
(8) 
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by time-reversal symmetry which generates a Z2 topological quantum number. 
In other words, the system (Eq. [6]) with M/2B > is a topologically nontrivial 
insulator which cannot be adiabatically tuned to a trivial insulator. More dis- 
cussions about the topological stability of this system will be presented in the 
next section. 

So far we have developed a simple 4-band lattice model for the Ei and Hi 
subbands and studied the properties of the bulk band structure and the low- 
energy critical point. We will now solve the model in a finite strip geometry to 
explicitly show the existence of topologically protected edge states for the non- 
trivial regime. Since the identity term e(fc) in the Hamiltonian (Eq. [S]) provides 
only an overall energy shift to all the states without changing the wavefunction 
profiles, we will omit this term below. The Hamiltonian (Eq.[6]) can be rewritten 
as 

n = Y^ {Asin(k^)T^ + Asixv{ky)T'^ + M{\<i)T'^) cj^Ck (10) 

k 

where F", a = 1 ... 5 form a Clifford algebra and Al(k) is defined in Eq. [6l For 
this model we have Fi = CT"=(8)s^ F^ = -ct^'^I, F^ = (t^(8)s^,F4 = a^'^sy, F^ = 
a^®l where ct' acts on the {Ei,Hi) space and acts on (-I-, — ) space. As stated 
above, around the F-point in the Brillouin zone this Hamiltonian behaves as a 
massive Dirac Hamiltonian in 2-d which means that for non-zero M the bulk 
spectrum will be gapped. Now we want to show that for certain values of the 
parameter M the system is in a topological insulator state characterized by a 
bulk energy gap and edge-states which lie in the gap. In order to study the 
edge-state spectrum we must first pick a particular edge. Since this simplified 
model has square lattice symmetry we pick the edges to be perpendicular to 
the j/-axis. This implies that ky is no longer a good quantum number since the 
system is not homogeneous in this direction, but instead has boundaries on the 
lines y = 0, L. We must Fourier transform ky into position space to obtain a 1-d 
tight-binding model which depends on the other momentum k^- We make the 
substitution 

j 

where k = k^. This gives the Hamiltonian: 

n = (12) 

M = Asai{k^)T^ -2B[2- MI2B -cos{k,^)]T^ (13) 
T = yF^ + BF^ (14) 

We expect edge states to be exponentially localized on the boundary so we 
look for solutions with the following ansatz [25,26] 

M) - A>„ (15) 
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where A is a complex number, j is the index of the j-th lattice site in the y- 
direction which runs from . . . L, and (f>a is a constant 4-component spinor with 
a — 1 . . A. Thus H^lja = Etpa gives 

A-lr„^0/3 + XT^p^[3 + Ma.0(pp = E^(3. (16) 

This Hamihonian is particle-hole symmetric, and time-reversal symmetric. Thus, 
at non-zero M a gap is opened around E — and at — we expect the edge 
states to have energy E — 0. So, at the F-point we generate the following simple 
equation 



y (A-i -X)T^ + B{X + A-i) + M(0)r^ 



= 0. (17) 



Multiplying both sides by gives 



^ (A-i - A) {^T^T^) <t>=- {B{X + X'^) + M{0)) ^ (18) 

The operator iF^F^ has eigenvalues ±1. First consider iT^T'^cj) = —cj), under 
which condition Eq. [18] becomes a quadratic equation for A which is easily 
solved to give 

-M{0)±^M^{0) + {A^-4B^) 
A(i,2) = . (19) 

It is trivial to see that if A is a solution for iT^T'^cf) = —0 then A^^ is a solution 
for the opposite eigenvalue iT^T^ip = (f). Denoting with s — 1,2 being 

the two eigenstates of iF''^F^ with eigenvalue -|-1(— 1), then a generic solution of 
the Schroedinger equation can be written as 



s s 

The open boundary condition can be expressed as ipa{j = 0) = 0. Since (/>s± 
are mutually orthogonal, this condition leads to as + bg — 0, Cs + dg — 0. 
On the other hand, the requirement of normalizability requires that only the 
exponentially decaying terms be included in ipU)- Consequently, a normalizable 
edge state solution exists only when |A(i)| < 1, |A(2)| < 1 (with Cg = ds = 0) or 
|A(i)| > 1, |A(2)| > 1 (with as=bs = 0). 

Such a derivation can be easily generalized to fc^^ 7^ case. Since [iF^F^, F-'^] = 
the term Asm{kx)T^ can be simultaneously diagonalized to add to the energy. 
The term involving cos(fea:) simply contributes to the mass term. It is conve- 
nient to choose the eigenstates 4>s± so that T^4's± — s0s±, in which condition 
the two edge state wavefunctions ips{j) and corresponding energy Es{kx) are 
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written as 



(-^(1) ~ ^(2)) '^'5+ 



(21) 



S 



Es{K) 



—As sm(kx) 



(22) 




-TO(fc, M) ± yfm{k,MY + ~4B2) 
-25(2 - M/2B - cos(fc^)). 



To have a normahzable solution wc still must have jA^^'^)] < 1 or lA^^'^)] > x, 
which leads to 



In other words, for a fixed M an edge state will exist only for the finite region 
of k that satisfies this inequality. 

For the dispersion in Eq. [22] we can see that two edge states with opposite 
eigenvalue are propagating in opposite directions. At k^ = 0, the edge states 
exist for < M < 45, that is, m(0, M) satisfies the inequality in Eq. [21] for 
M/2B e (0, 2). When k^ is varied from to tt, the edge states are dispersive 
until the inequality (Eq. I23|) is violated at some finite kx^ which, physically, 
is the wavevector where the edge states touch the bulk (extended) states. In 
Fig. [3] we have plotted the energy spectrum as a function of k^ from (numerical) 
exact diagonalization, and from our analytic solution. In our analytic solution, 
and with fixed M, we solve for the A:a;-points where the edge state solution is no 
longer normahzable, which occurs when the edge state becomes degenerate with 
the bulk states. It is at these fc-points where the edge state dispersion enters 
the bulk and is no longer visible. 

4 Topological Properties of the Helical Liquid 

Similar to the case of quantum Hall effect, the topological properties of the 
QSH system can be understood by studying the low energy effective theory of 
the edge states. The edge states of the quantum Hall effect are described by a 
chiral Luttinger liquid theory [27]. Such a chiral theory only contains, say, left 
movers, so that the current carried by the edge excitations cannot be destroyed 
by any perturbation due to the absence of counter-propagating modes in which 
to backscatter. In comparison, the effective theory of the quantum spin Hall 
edge states do contain both left and right movers, and can be written as 



where -ffpcrt represents possible perturbation terms, including disorder and/or 
electron-electron interactions. 

If time- reversal symmetry is not present, a simple "mass" term can be added 
in iJpert so that the spectrum becomes gapped: 



2B < m{k, M) < 2B 



(23) 




(24) 
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Figure 3: (Upper)Exact diagonalization of 1-d tight binding model. (Lower) 
analytic solution for bulk (blue) and edge states (red). Bulk bands are plotted 
vs kx for several different values of ky in order to reproduce the look of the exact 
diagonalization results. 



However, the time-reversal symmetry of the electron system is expressed as 

T-Vfe+T = V-fc,-, T-Vfc-T = (25) 

which implies 

Consequently, i?mass is a time-reversal symmetry breaking perturbation. More 
generally, if we define the "chirality" operator 



C = iV I - 7V_ = 



dk 

2^ 



any operator that changes C by 2(2n — 1), n e Z is time- reversal odd. In other 
words, time-reversal symmetry requires that the Hamiltonian i?port only include 
processes of 2n particles back scattering, such as ^ipp-ipp' - . Therefore, 

the most relevant perturbation in normal metals, i.e., ipk+'^'k'-, is forbidden by 
time- reversal symmetry. This is essential for the topological stability of the edge 
states. This edge state effective theory is thus non-chiral and is qualitatively 
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different from the usual spinless or spinful Luttinger liquid theories. It can be 
considered as aa example of a new class of one-dimensional critical theories, 
describing a "helical liquid" [13,14]. Specifically, for the non-interacting case 
no time-reversal invariant perturbation is available to induce back scattering, 
so that the edge state is robust under any disorder. 




Figure 4: (a) The energy dispersion of a one-dimensional time-reversal invariant 
system. Kramer's degeneracy is required at fc = and A: = tt, so that the energy 
spectrum always crosses the Fermi level 4n times, (b) The energy dispersion 
of the helical edge states on one boundary of the QSH system (solid lines). At 
fc = the edge states are Kramers' partners of each other, while at fc = tt they 
merge to the bulk and pair with the edge states of the other boundary (dash 
lines). In both (a) and (b), red and blue lines form Kramers' partners of each 
other. 

Now consider the case with two flavors of helical edge states on the bound- 
ary, i.e., a 1-d theory with two left movers and two right movers, with the 
Hamiltonian 

s=l,2 

Due to the same reason as the one-flavor case, the single particle back scattering 
terms between Kramers' doublets V'L-i-''/'-fc s- forbidden. However, the back 
scattering between different flavors such as V'li+V'-fc 2- '^^^ tie non-vanishing, 
which makes such a theory qualitatively equivalent from a spinful Luttinger 
liquid with four degrees of freedom. A mass term like / (^4'ii+''Pk2- + h.c.^ 
can immediately open a gap in the system, which implies that the existence of 
gapless edge states is not generic, as shown in Fig. [S] In other words, two copies 
of the helical liquid becomes a topologically trivial theory. More generally, an 
edge system with time-reversal symmetry is a nontrivial helical liquid when 
there is an odd number of left (or right) movers, and trivial when there is an 
even number. Such a property makes it natural to characterize the topology in 
the QSH systems by a Z2 topological quantum number. 
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Figure 5: Brown areas are regions of bulk states. Red and blue dispersing 
edge states represent Kramers' partners (a)An even number of pairs of fermion 
branches crossing at k points which do not transform to themselves under time- 
reversal (b) A slight perturbation added to the system causes the edge states 
in (a) to form a gap. Kramers' theorem is still satisfied but these edge states 
are unstable to gap formation (c) Two pairs of fermion branches that initially 
crossed at the special point k = are shown after an infinitesimal perturbation is 
added. A gap is formed and Kramers' theorem is satisfied but this configuration 
is also unstable to gap formation (d)A single pair of fermion branches crosses at 
k = 0. A perturbation cannot open a gap because in that case there would be 
two states which were singly degenerate which will not satisfy Kramers' theorem, 
thus this configuration is stable. 



There is also an alternative way to understand the qualitative difference 
between even and odd branches, which is expressed as a "no-go" theorem in 
Ref. 13: There is always an even number of Kramers' pairs at the Fermi en- 
ergy for an arbitrary one-dimensional band structure. Such a no-go theorem 
is a time-reversal invariant analog of the Nielsen-Ninomiya no-go theorem for 
chiral fermions on a lattice [28,29]. For spinless fermions there are always an 
equal number of left movers and right movers at the Fermi level, which leads 
to the fermion doubling problem in odd spatial dimensions. Similarly, for a 
time-reversal symmetric system with half-odd integer spin, Kramers' theorem 
requires each energy band to be accompanied by its Kramers' partner, so that 
the number of low energy channels are doubled again. As shown in Fig. 2] (a) , 
a Kramers' pair of states at fc = must re-form into pairs with each other when 
k goes from to tt and 27r, which requires the bands to cross the Fermi level An 
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times. Seemingly disheartening, there is actually an exception to this theorem 
analogous to the reason why a chiral liquid exists in the quantum Hall effect. 
The exception is that a helical liquid with an odd number of fermion branches 
can occur if it is holographic, i.e., appears at the boundary (edge) of a 2d bulk 
system. As shown in Fig. |31 in this case the edge states are Kramers' partners 
of each other at fc = 0, but they merge into bulk states at some finite kc, so 
that they don't have to be paired with each other at = tt. More accurately, at 
k > kc the edge states on both the left and right boundaries become bulk states 
at fc > fee, and form a Kramers' pair with each other, as shown in Fig. 3] (b). 

The no-go theorem also provides a physical understanding of the topological 
stability of the helical liquid. Any local perturbation on the boundary of a 
2d QSH system can be considered as coupling a "dirty surface layer" with the 
unperturbed helical edge states. Whatever perturbation is considered, the "dirty 
surface layer" is always a one-dimensional system, so that there are always an 
even number of Kramers' pairs of low energy channels. Since the helical liquid 
has only an odd number of Kramers' pairs, the coupling between them can only 
annihilate an even number of Kramers' pairs, so that at least one-pair of gapless 
edge states can survive. 

Similar to the correspondence of chiral edge states and bulk Chern number in 
quantum Hall effect [9] , the Z2 topological stability of helical edge states can be 
determined by a Z2 topological quantum number defined in the two-dimensional 
Brillouin zone, as shown in Refs. 15 and 30. Several alternative ways to define 
the Z2 invariant and its generalization to three-dimensional topological insu- 
lators were also proposed in the literature [31-36]. In the quantum Hall case, 
the Chern number can be defined by introducing twisted boundary conditions, 
even when disorder and interactions are present so that momentum is not good 
quantum number [37]. Similar twisted boundary conditions can be introduced 
to define the Z2 invariant when disorder is considered [38]. However, the gen- 
eralization to the many-body case with interactions is not so straightforward, 
since the Kramers' degeneracy is only defined for a single electron state. Some 
efforts have been made toward this direction recently [39]. 

In summary, in this section we have discussed the topological stability of the 
edge states, which is essential for the quantum spin Hall effect in HgTe quan- 
tum wells to be realized in an experimental system. Although the edge states 
are non-chiral, one can still detect them through the longitudinal conductance 
it contributes. In the ballistic limit, a helical liquid provides a metallic one- 
dimensional channel, which leads to a conductance of e^//i. For a system with 
two boundaries the net longitudinal conductance will be 2e^//i. In the rest of the 
paper we will review the experimental observations of the transport properties 
of HgTe quantum wells and present theoretical analysis of the results. 

5 Quantum Well Fabrication 

As discussed above, we anticipate that the transition from a normal to an in- 
verted band structure in HgTe QWs coincides with a phase transition from a 
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trivial insulator to the QSH insulator. Experimentally, we have investigated the 
transport properties of HgTe samples with QW width in the range from 4.5 nm 
to 12.0 nm, so as to cover both the normal and the inverted band structure 
regime. The samples were grown by molecular beam epitaxy (MBE) and the 
layer sequence is schematically depicted in Fig. [S] The wells are sandwiched 
by Hgo.3Cdo.7Te barriers and n-type modulation doped using I-doping on both 
sides of the QW layer. Recent advances in the growth of HgTe-based QW struc- 
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Figure 6: Schematics of the layer sequence of the MBE-grown quantum well 
structures. 

tures, that were incorporated in these structures are documented in Ref. [40]. 
For example, increasing the spacer width between the quantum well and the 
doping layer results in an enhancement of the carrier mobility /i, and samples 
with mobilities of several 10^ cm^/(Vs) even at low densities n < 5 x 10^-'^ cm~^ 
have been used for the actual measurements. In such samples, the elastic mean 
free path is of the order of several microns. 

The devices have been structured by means of optical and electron beam (e- 
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beam) lithography. So as to avoid Hg evaporation, the HgTe layer has to be grow 
at temperatures as low as 180°C and the sample has to be kept well below this 
temperature during the entire process of fabrication. Standard processes based 
on PMMA (polymethylmethacrylate) as a resist can not be used, since these 
usually demand temperatures up to 200° C for bake-out. While we fabricated 
our first HgTe-based nanostructures using a standard photoresist for the e-beam 
steps [41], we have, for the samples discussed here, employed a dedicated low 
temperature PMMA technology. In combination with Ar-ion etching we can 
reproducibly fabricate structures with dimensions down to ca. 100 nm. 

For the investigation of the QSH effect, devices in a Hall bar geometry of 
various dimensions have been fabricated from QW structures with a QW width 
of 4.5 nm, 5.5 nm, 6.4 nm, 6.5 nm, 7.2 nm, 7.3 nm, 8.0 nm and 12.0 nm, 
respectively. The layout of the devices is sketched in Fig. [T] 



I 
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Figure 7: This schematic sketch of the devices shows the QSH edge states under 
the gate (shaded region). The ungated parts of the sample are n-type. 

A 110 nm thick layer consisting of a Si3N4/Si02 superlattice is deposited on 
top of the structure by plasma-enhanced chemical vapor deposition (PECVD). 
This layer serves as a gate insulator between the semiconductor and the (Ti/Au) 
gate electrode on top of the structure. The top gate covers the leads of the Hall 
bar only to a small extent, while the largest parts of the leads are not gated, 
which implies that the mesa is always contacted by n-type HgTe leads. The 
gate insulator has excellent break-through properties and allows gate voltages 
in the range from -5 V to -1-5 V which is sufhcicnt to induce a large change in 
carrier concentration [40,42]. However, in most structures not the full voltage 
swing can be applied. Hinz et al. have reported hysteresis effects in HgTe QWs 
at extreme gate voltages [43]. These effects are attributed to the filling and 
emptying of trap levels between the insulator and the semiconductor. In the 
experiments described below, care was taken to keep the gate voltages in the 
non-hysteretic regime. 
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6 High Field Characterization 



For the investigation of the QSH effect, samples with a low intrinsic density, 
i.e., n{Vg = 0) < 5 X 10^^ cm"^, are studied. Transport measurements are done 
in a ■^He/^He-dilution refrigerator with a base temperature T < 30 mK and 
a *He cryostat {T — lA K) fitted with a vector magnet system, which allows 
for magnetic fields up to 300 mT in variable direction. When a negative gate 
voltage Vg is applied to the top gate electrode of the device, the usual decrease 
of the electron density is observed. In Fig. [8] A, measurements of the Hall 
resistance R^y are presented for a Hall bar with length L = 600 fim and width 
W = 200 /xm, respectively. The decrease of the carrier density is reflected in an 
increase of the Hall coefficient, when the gate voltage is lowered from V. For 
voltages down to -1.3 V, the density is decreased linearly from 3.5 x 10^"'^ cm^ 
to 0.5 X 10^^ cm^ (Fig-ElB). For even lower gate voltages, the sample becomes 




Figure 8: A: Hall resistance Rxy as measured for various gate voltages, indicating 
the transition from ri- to p-conductance. B: The gate-voltage dependent carrier 
density deduced from the Hall measurements. 

an insulator, with the Fermi level in the band gap. Remarkably, when a large 
negative voltage Vg < —2 V is applied, the sample becomes conducting again. 
It can be inferred from the sign change of the Hall coefficient that the device 
is now p-conducting. Obviously, the Fermi level has now been pushed through 
the entire gap and into the valence band. 

The transition from n- to p-conductance is also reflected in the behavior of 
the longitudinal resistance Rxx as a function of Vg (cf. Fig. [T3|) . For a decreasing 
electron density the resistance rises until it reaches a maximum when the Fermi 
level is in the gap. When the Fermi energy finally is pushed into the valence 



17 



band and the sample is p-conducting, R^^ decreases again by some orders of 
magnitude. A similar transition from n- to p-conductance via an intermediate 
high resistance regime has been observed for a variety of QW structures with 
low intrinsic densities. The exact value of the maximum resistance, however, 
depends critically on the nature of the band structure, and thus on the actual 
Q W width of the individual devices (cf. Section 7, below) . 

A similar signature in transport for a transition from n- to p-conductance is 
possible and has been reported for the zero-gap material graphene [44] . However, 
since this material is a semi-metal, no insulating regime is observed between n- 
and p-conductance. Thus, the QSH effect is experimentally not accessible in 
graphene even though it was explicitly predicted for this material [11]. It was 
later understood theoretically that a possible energy gap in graphene due to 
spin-orbit coupling is far too small to be experimentally observable [16, 17]. 

The peculiar band structure of inverted HgTe quantum wells gives rise to a 
unique Landau level (LL) dispersion. For a normal band structure, i.e., dqw < 
dc, all Landau levels are shifted to higher energies for increasing magnetic field 
[Fig. [5] (a)]. This is the usual behavior and can be commonly observed for a 
variety of materials. When the band structure of the HgTe QW is inverted 
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Figure 9: Landau level dispersion for (a) a 4.0 nm QW and (b) a 15.0 nm 
QW, respectively. The qualitative behavior is representative for samples with a 
normal and an inverted band structure, respectively. 

{dgw > dc), the LL dispersion is markedly different [Fig. ^ (b)]. Due to the 
mixing of the electron-like and hole-like states, one of the states of the Hi 
subband is a pure heavy-hole state (M = —3/2). Consequently, the energy 
of the corresponding Landau level decreases with increasing magnetic field. In 
adition, one of the valence band LLs has mainly electron character and, thus. 
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shifts to higher energies with magnetic field. This leads to a crossing of these 
two peculiar LLs for a finite magnetic field. The exact value of the crossing field 
-Across depends on the QW width. Fig. [TUl shows theoretically calculated S- field 
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Figure 10: The S-field value of the LL crossing determined experimentally (blue 
squares) and calculated theoretically (red circles). For the energy gap Eg (black 
triangles), the negative values indicate the inverted band structure. 

values for the Landau level crossing for quantum wells with an inverted band 
structure as red circles. The existence of such a LL crossing is a clear indication 
for the occurrence of an inverted band structure; for narrow QWs (dgw < dc) 
with a normal band structure, there is no LL crossing. The figure also plots 
the energy gap of the HgTe QWs, defined as Ei — Hi, as black triangles. In 
narrow wells (dgw < dc) with a normal band structure, our definition yields 
a positive gap, whilst for inverted wells we have an increasingly negative gap. 
For dQw ~ 8.5 nm, the Ei subband falls below the H2 band (cf. Fig. (2), and 
as a result the absolute value of the gap decreases subsequently for larger QW 
widths. 

The crossing of the conduction and valence band Landau levels can be ver- 
ified experimentally by quantum Hall experiments. Fig. [TT] (a) shows quantum 
Hall data for an inverted quantum well with a width of 6.5 nm. For gate volt- 
ages Vg > —1.0 V and Vg < —2.0 V, the slope of the Hall signal directly shows 
that the Fermi level is firmly in the conduction and valence band, respectively. 
When the Fermi level is adjusted in the gap, i.e., —1.9 V< Vg < —1.4 V, the 
Hall resistance shows insulating behavior at zero field. 

However, at higher magnetic fields the Hall resistance for these gate voltages 
exhibits a re-entrance of the n = 1 quantum Hall plateau. This is a direct 
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Figure 11: (a) Hall resistance, R^y, of a. {L x W) = (600 x 200) A^m^ QW 
structure with 6.5 nm well width for different carrier concentrations obtained 
for different gate voltages Vg in the range of -1 V to - 2 V. For decreasing Vg, the 
n-type carrier concentration decreases and a transition to a p-type conductor is 
observed, passing through an insulating regime between -1.4 V and -1.9 V at 
B = T. (h) The Landau level fan chart of a 6.5 nm quantum well obtained 
from an eight band k-p calculation. Black dashed lines indicate the energetic 
position of the Fermi energy, Ep, for Vg = —1.0 V and —2.0 V. Red and green 
dashed lines correspond to position of the Fermi energies of the red and green 
Hall resistance traces of (a). The crossing points of the Fermi level with the 
respective Landau levels are marked by arrows of the same color. 
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consequence of the crossing of the Fermi level with the conduction-band derived 
LL for the green trace {Vg = —1.4 V) and with the valence band LL (e.g, for the 
red trace where Vg — —1.8 V). As shown Fig. [TT] (b), the experimental results 
are in good agreement with the theoretically calculated LL dispersion. The 
crossing point of the Landau levels in magnetic field, Across, can be determined 
accurately by tuning the Fermi level through the energy gap. This method has 
been used for various samples with an inverted band structure; the results are 
included in Fig. [TO] (blue squares) . 

The observation of a re-entrant quantum Hall state is a clear indication of 
the non-trivial insulating behavior, which is a prerequisite for the existence of 
the QSH state. In contrast, trivial insulating behavior is obtained for devices 
with dgw < dc- For a normal band structure, the energy gap between the 
lowest Landau levels of the conduction band and the valence band, respectively, 
increases in magnetic field [cf. Fig. [9] (a)] . Thus, a sample remains insulating in 
magnetic field, if the Fermi energy is located in the gap for B ^ 0. 

The physics of this re-entrant quantum Hall state can basically be under- 
stood within the simple 4-band model. We will consider our model given in 
Eq. [6] expanded near k = and in the presence of a uniform i3-field {eB > 0) in 
the z-direction. We choose the symmetric gauge with A = ^{—y,x). We make 
the substitution k ^ (k -I- eA) and define new operators 

»+ = ''('•'+ + ^«) w 

where k± = [k^ ± iky) and z = x + iy. These operators obey the commutation 
relations [7r+, 7r_] = -jr^ with the magnetic length £b ~ (h/eB)^^'^. Using these 
commutation relations we can define raising and lowering operators 

a^-n^ , = (28) 

[a, at] = 1. (29) 
Using these operators we rewrite our Hamiltonian as 

«=('^'o'°\.(!.„)) w 



I tg 



V2A, . _ 

— - — (a' a +aa'^) 



I iB 
^ z 

+ 2^^ 



(31) 
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where cr^ = |(cr^ ± ia"^). It should be noted in this context that the material- 
specific parameter B of the Dirac model, which plays a crucial role in describing 
the re-entrant quantum Hall effect, is not related to the magnetic field. The spec- 
trum of this Hamiltonian can be solved since only a finite number of harmonic 
oscillator Landau levels are coupled. We find the energy spectrum 



B 



= -uj^n -a-^±J —2— + (M - nw^ - -Lolf)^ (32) 



2 ^ V 



where ujq = -J-i?, oj^^ — -^D, a — ± and n = 0, 1, 2, . . . This spectrum has 
"zero modes" given by 

K = -M-hu:^ -u^^) (33) 



2^ 

2' 



= M+hu^^ +u:^) (34) 



for M/2B < and 



El = M-\{ujE+oj^) (35) 
E'_ = -M+\{u^E+^^) (36) 

for M/2B > 0. For M/2B < one can see that the set of zero-modes never 
cross. However, for M/2B > 0, i.e., in the inverted QSH regime, the zero 
modes do cross at a field of Across — f "bTd' ^'^^ typical parameters in the 
inverted regime [18, 45] the critical field is just a few Tesla. Since the Dirac 
mass M reflects the energy gap, the linear dependence of Across on M in the 
4-band Dirac model is in good agreement with the results obtained by the 8x8 
Hartree calculations shown in Fig. [TOl The typical Landau level spectra for the 
non-inverted and inverted regime, respectively, presented in Fig. [12] (a), (b) are 
obtained by the simple 4-band model. 

As illustrated in Fig. [T^] there are four different regimes of interest. Two have 
quantized Hall conductances of zte^/h and the other two have vanishing Hall 
conductances. Case (i) is described hy B > Be, while the chemical potential 
is placed between the zero modes. Here no edge states exist and this yields 
<7xy = 0. Identical behavior is obtained for a trivial insulator independent of 
the S-field strength, if the chemical potential is in the energy gap. Case (ii) also 
shows axy = 0. However, there should be counter-propagating edge states that 
carry opposite Hall current, giving rise to the QSH effect. When the chemical 
potential is slightly below the two zero modes as in case (iii), only one of them 
will provide edge states. Thus, there is a single state at each edge, resembling 
the well-known quantum Hall regime with axy = e^/h. A similar result is 
obtained for case (iv), where the chemical potential is just below the zero modes. 
However, since the single edge state now originates from the other zero mode, 
the Hall conductance changes sign, i.e., a^y = —e^/h. If the chemical potential 
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is fixed and unequal to the energy of tfie LL crossing i? (Across), an increase of 
the magnetic field gives rise to the re-entrant quantum Hall state. While the 
Hall conductance vanishes for case (i) and (ii), an intermediate regime with 
axy = ±e^//i according to case (iii) and (iv), respectively, will be obtained, 
while the sign of a^y depends on the exact chemical potential. 

7 Experimental observation of the QSH effect 
at ^ = 

The existence of the QSH effect is revealed when small Hall bars {L x W = 
20.0 fim X 13.3 ^m) are studied. For QW devices with a normal band structure, 
the sample shows trivial insulating behavior (Fig. [13]). The resistance of several 
M^l for the insulating regime can be attributed to the intrinsic noise level of 
the measurement setup, so that the conductance basically vanishes. For a de- 
vice with an inverted band structure, however, the resistance stays finite, not 
exceeding 100 kfi. This behavior is reproduced for various Hall bars with a QW 
width in the range from 4.5 nm to 12.0 nm: while devices with a normal band 
structure, i.e., dqw < dc ~ 6.3 nm, show trivial insulating behavior, a finite 
conductance in the insulating regime is observed for samples with an inverted 
band structure. We now proceed ro relate this observed finite conductance to 
the predicted edge state behavior. In the theoretical modelling of Ref. 18, ideal 
contacts to the QSH edge states were assumed, i.e., the potential of the edge 
states can be probed without any influence due to the contacts. However, in 
reality the leads are always n-type for the investigated samples and remain such 
independently of the applied gate voltage. Unlike for a quantum Hall system, 
the edge states in a QSH system propagate in both directions at a given edge 
(cf. Fig. [7]). Thus, the edge states entering a contact carry a different potential 
due to the different sources. When the edge states enter an n-type contact, they 
necessarily equilibrate with the bulk states and backscattering becomes possible. 
Hence, the voltage difference between neighboring contacts does not vanish even 
though they are connected by non-dissipative edge channels. In other words, 
a voltage probe induces additional resistance, even when no voltage is applied 
to it. Such a dissipation occurs due to the large size of the n-type contact, 
which implies a huge number of channels the incoming electron can scatter to. 
Consequently, an infinitesimal dissipation in the contact can cause complete de- 
coherence between the incoming and outgoing electrons. This analysis can be 
easily confirmed by a Green's function calculation of the ballistic conductance in 
the Keldysh formulism. In contrast, in a quantum Hall system, a voltage probe 
does not induce additional dissipation, because all the states are propagating to 
the same direction. The de-coherence effect of a contact, although also it exists, 
does not affect the current through it. 

A simple Landauer-Biittiker formalism [46] can be applied to determine the 
influence of the n-type contacts. The current Ii in a contact i (cf. Fig. [7|) can 
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Figure 12: (a) Landau level fan diagram for the El and HI subbands in the non- 
inverted regime. (b)Landau level fan diagram for the El and HI subbands in the 
inverted regime. Notice that the zero-modes from each set of subband Landau 
levels cross at a finite yB-field. Lower case roman numerals indicate separate 
regions in energy and magnetic field. (c)-(f)Schematic energy spectrum of the 
zero-modes and edge states for region i,ii,iii, and iv respectively. (c)Region 
i produces a vanishing hall-conductance, (d) Region ii produces a vanishing 
hall-conductance but the Fermi-level passes through two branches of quantum 
Hall edge states which are oppositely circulating and carry opposite current, 
(e) Region iii produces a Hall-conductance of +e'^/h, and the Fermi-level only 
crosses one of the sets of edge states, (f ) Region iv produces a Hall-conductance 
of —e^/h, and the Fermi-level only crosses one of the sets of edge states. 
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Figure 13: The longitudinal resistance of a 4.5 nm QW (a) and an 8.0 nm QW 
(b) as a function of gate voltage, respectively. 

be calculated from 

It = ie/h)^T,j{^ij - fi,). 

3 

Since the edge states are non-interacting, the transmission coefficient Ty is unity 
for neighboring contacts and zero otherwise. In the measurements, a current I 
is injected from contact 1 to contact 4, i.e., Ii ~ —I4 — I, while = vanishes 
for all other contacts serving as voltage probes. When total equilibration of all 
edge states in each contact is assumed, a four-terminal resistance for neighboring 
voltage probes of (/i/2e^) is obtained. The two-terminal resistance is determined 
by the number n of the voltage probes between the current source and drain. 
Between each pair of neighboring contacts, a voltage Vi = I-{h/e^) drops, adding 
up to a total voltage drop of (n -I- l)Vi. Thus, the two-terminal resistance is 
(3/i)/(2e^) for a Hall bar geometry with two voltage probes on each side as 
shown in Fig. \7\ 

However, the obtained finite resistance due to the n-type contacts still does 
not give a quantitative explanation for the observed resistance R ~ 100 kfi, 
which is significantly higher than four terminal resistance {h/2e?) « 12.9 kfJ 
one anticipates for the geometry used in the experiments. The enhanced resis- 
tance in these samples with a length of L = 20 /im can be readily understood 
when inelastic scattering is considered. While, as discussed above, the heli- 
cal edge states are robust against single-particle elastic backscattering, inelastic 
mechanisms can cause backscattering. 
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For n-doped HgTe quantum wells, the typical mobility of the order 10^ cm^/Vs 
implies an elastic mean free path of the order of 1 /zm [41]. Lower mobilities 
can be anticipated for the QSH regime. The inelastic mean free path, which 
determines the length scale of undisturbed transport by the QSH edge states, 
can be estimated to be several times larger due to the suppression of phonons 
and the reduced electron-electron interaction at low temperatures. Thus, the 
inelastic scattering length is of the order of a few microns. 

For the observation of the QSH conductance, the sample dimensions were 
reduced below the estimated inelastic mean free path. When Hall bars with 
a length L — 1 /im are studied, a four-terminal resistance close to h/{2e^) is 
observed. The threshold voltage Vth is defined such that the QSH regime is in 
the vicinity of = Vth- The slight deviation of R from the quantized value 
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Figure 14: The longitudinal resistance is measured as a function of the gate 
voltage for two devices with L = 1/im. The width is 1 /xm (black and blue) 
and 0.5 /xm (red), respectively. The black and red trace were obtained at 1.8 K, 
the blue one at 4.2 K. 

/i/(2e^) can be attributed to some residual scattering. This is an indication 
that the length of the edge states still exceeds the inelastic mean free path. 
However, a significant increase of the conductance is obtained by reducing the 
sample size, which demonstrates clearly the effect of scattering on the transport 
by QSH states. 

The results presented in Fig. [H] provide evidence that transport in the QSH 
regime indeed occurs due to edge states. The two devices with W = 1.0 /im 
and W — 0.5 /im, respectively, were fabricated from the same QW structure. 
The resistance values of the two devices differ significantly for the n-conducting 
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regime, where transport is determined by bulk properties. For the insulating 
regime, however, both devices show the same resistance, even though the width 
of the devices differs by a factor 2. For all three measurements, no full transition 
to the p-conducting regime can be observed. This is attributed to charging of the 
quantum well layer during the e-beam process required to fabricate such small 
samples. The fluctuations, which are observable in all traces, are reproducible. 
Thus, they must have a physical reason, which is not fully clear yet, and can 
not be caused, e.g., by electrical noise. The pure QSH signal, however, is robust 
against an increase of temperature as long as /csT is much smaller than the 
two-dimensional bulk energy gap. 



8 Suppression of the quantum spin Hall effect 
by a magnetic field 

Further clear evidence for the QSH effect is obtained by measurements in a 
finite magnetic field, which destroys the time reversal symmetry of the QSH 
states. In order to observe this effect, the dependence of the QSH conductance 
is studied in a magnetic field in various directions with respect to the plane 
of the 2DEG. For this purpose an inverted QW in Hall bar geometry with 
(L X W) = (20.0 X 13.3) firn^ was investigated at 1.4 K in a vector magnet 
system. In a perpendicular field configuration the QSH conductance decreases 
rapidly when the magnetic field is applied. A cusp-like conductance peak is 
observed with a full width-half maximum field Sfwhm of 28 mT (Fig. fTS]) . 
which is even smaller at lower temperatures. At T = 30 mK a typical width of 
'Bfwhm = 10 mT is observed. When the magnetic field direction is rotated into 
the plane of the 2DEG the conductance peak widens strongly, reaching a value 
of iJpwHM ~ 0.7 T for the parallel configuration. 

From the results shown in Fig. [15] it is evident that a perpendicular field has 
a much larger influence on the helical liquid than an in-plane field. 

Such a large anisotropy can be understood by a slightly modified version of 
the four-band effective model we discussed in Sec. 3. To first order in pertur- 
bation theory, the gap induced by a perpendicular and an in-plane field can be 
determined by the matrix elements of the corresponding operator: 

Sgap± = |(fc = 0,+| (z-r X J+/iBA4)s^|/^ = 0,- 
Sgap|| = \{k = 0,+\^iB{M.B^ + MyBy)\k = 0,-)\ (37) 

where r, j are electron position and current operators, respectively, and z is 
the unit- vector perpendicular to the 2DEG plane. The states j/c = 0, ±) are the 
two edge states propagating toward opposite directions on the same boundary. 
The Zeeman coupling matrix elements M^- y ^ in the four-band effective model 
can be determined by standard perturbation procedures from the original Kane 
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Figure 15: Four-terminal magnetoconductance, Gi4,23i in the QSH regime as 
a function of tilt angle between the plane of the 2DEG and applied magnetic 
field for a d = 7.3 nm QW structure with dimensions {L x W = 20 x 13.3) /xm^ 
measured in a vector field cryostat at 1.4 K. 



model, which has the form [45] 



M,, = 



( 9EA 



\ 



'9E1^ 



-9H± I 

f n \ ( 

,My = 



V 



.911 



\ 



i9\\ \ 
J 



(38) 



Since the edge state wavefunctions of the effective model have been obtained 
in Sec. 3, the gaps -Bgapi and -Egap|| can be calculated for given g-factors and 
Zeeman couplings Sx.y,z- However, there is a subtlety that has to be considered 
before carrying out this calculation. The effective model (Eq. [1]) consists of 
two decoupled blocks, which means the coupling between |i51, , \H1, and 
|i?l,— ), |iJl,— ) is ignored. Thus for such an effective model, the current 
operator j = Vki?(k) is also block diagonal. Since the wavefunction of the edge 
state \k, )) consists only \El, +{—)) and \H1, -|-(— )) components, the orbital 
magnetization operator z • f x j has a vanishing matrix element between them, 
which means the orbital effect of the magnetic field does not contribute to the 
gap -Egap^. However, this turns out to be an artifact of the simplified model, 
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which comes from the absence of bulk inversion asymmetry (BIA) terms. When 
the breaking of bulk inversion symmetry is taken into account, an additional 
term is induced in the four-band effective model [45] : 



Hbia 



I 



\ -A 



A 



-A \ 



(39) 



which describes the mixing between )) and (+)). Such a term 

preserves time-reversal symmetry, and does not close the bulk gap, and thus 
does not affect the topological stability of the nontrivial insulator phase. Con- 
sequently, all the theoretical discussions based on the four-band model that we 
presented in the earlier part of the paper are only slightly modified. Never- 
theless, such a term becomes important when we consider the effect of per- 
pendicular magnetic fields. Since the bulk- inversion symmetric terms have 
a vanishing contribution to the edge state gap, the contribution from BIA 
terms becomes the leading order. The BIA term A and all other parame- 
ters A, B,C, D, M, gE±, gH±, g\\ in the model can be obtained from envelope 
function calculations of the quantum well [45] . For the quantum well thickness 
d — 7.0 nm, these parameters are given as 



A{eY ■ A) 


B(eV- A^) 


^(eV) 


M(eV) 


A(eV) 


9E± 


9H± 


.911 


3.645 


-68.6 


-51.2 


-0.010 


0.0016 


22.7 


-1.21 


-20.5 



By inputing these parameters into the tight-binding model (Eq. ^ the edge 
state wavefunction can be obtained numerically. For a magnetic field of IT, 
the gap induced is i?gap_L = 3.1 meV if it is perpendicular to the 2DEG plane, 
or fi'gapii — 0.3 meV if parallel. Since the perpendicular and in-plane Zeeman 
coupling are of the same order, we see that the main contribution of i?gap± 
comes from the orbital effect of the magnetic field. Compared to the Zee- 
man contribution, the orbital effect is equivalent to an "effective g-factor" of 
order 10^. Besides the numerical calculations, the reason for such a large or- 
bital effect can also be understood by a rough order-of-magnitude estimation. 
Since the current operator still doesn't mix \E{H), +) and \E{H), — ) even with 
BIA terms in the Hamiltonian, the only contribution to -Bgapi comes from 
the weight of \E{H), — ) components in the edge state wavefunction \k = 0, 
which is proportional to A/|M| from first order perturbation theory. Thus 
we can estimate the matrix element of the orbital magnetization operator as 
(fc = 0, -|-| z • f X j |fc = 0, — ) ~ {A/\AI\) evF^, where we use evp and the edge 
state width ^ as the estimates of operators j and f, respectively, vf and ^ of the 
edge states can be estimated by = A/h, ^ ~ ^/|M|. Thus £;orb - f 
Compared with the Bohr magnon we have the "effective ^-factor" as 



e AA2 



1 



2'meVp 



h JVP /is Af2/A 



(40) 
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Since vp = A/h ~ 5.5 x lO^m/s, wc obtain gcs is around 50. Such an analysis 
also leads to the prediction that the anisotropy between perpendicular and in- 
plane magneto-conductance is further enhanced when the thickness goes closer 
to the critical thickness dc from the inverted side, due to the decrease of \M\. 

9 Conclusion 

Wc have reviewed our current theoretical understanding of the QSH state, with a 
focus on the experimental realization in HgTe quantum wells. We discussed the 
electronic structure of the HgTe quantum wells in terms of a simple 4 band model 
which contains the essential physics. Wc demonstrated the topological quantum 
phase transition as the quantum well thickness cIqw is varied and we explicitly 
showed the analytic solution of the helical edge states on the topologically non- 
trivial side of the phase transition. Wc discussed the topological stability of the 
helical edge states using the concrete example of the 4 band model for HgTe 
quantum wells and presented the experimental realization of the QSH effect 
in HgTe QW structures for various QW widths. These experiments clearly 
demonstrated the edge channel character of the QSH effect and that the inverted 
band structure, which occurs for HgTe QW with cIqw > 6.3 nm, is essential. 
Deviation from the quantized conduction value for large samples and the effect 
of an applied magnetic field could be related to back-scattering introduced by 
time-reversal symmetry breaking processes. 

The actual experimental realization of the QSH effect opens up the opportu- 
nity for investigations of new theoretical concepts [47] , concerning fundamental 
aspects and applications, utilizing the spin polarized properties of the helical 
edge channels. 
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and by the US Department of Energy, Office of Basic Energy Sciences under 
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